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We calculate the critical temperature and free energy of the gluon plasma 
using the dilaton potential [1] in the gravity theory of AdS/QCD. The 
finite temperature observables are calculated in two ways: first, from the 
Page-Hawking computation of the free energy, and secondly using the 
Bekenstein-Hawking proportionality of the entropy with the area of the 
horizon. Renormalization is well defined, because the T = theory has 
asymptotic freedom. We further investigate the change of the critical 
temperature with the number of flavours induced by the change of the 
running coupling constant in the quenched theory. The finite temperature 
behaviour of the speed of sound, spatial string tension and vacuum expec- 
tation value of the Polyakov loop follow from the corresponding string 
theory in AdS 5 . 

1 INTRODUCTION 

In the bottom-up approach of AdS/QCD important properties of pure glue QCD are 
encoded in a phenomenological gravity theory through the introduction of a dilaton po- 
tential. The fifth dimension plays the role of an inverse energy scale, which necessitates 
that the dilaton is not constant, but runs with this scale. In a previous paper [1] we have 
fixed the ultraviolet behaviour of the potential to the two loop beta function of QCD and 
parametrized the infrared part in such a way that the heavy QQ potential is reproduced. 
In five dimensions the string connecting the Q and Q hangs into the bulk fifth dimension 
and thereby is sensitive to the geometry of the five dimensional space. An obvious next 
step is to consider properties of QCD in other simple environment where the geometry of 
the five dimensional space changes in a controlled manner. 

Finite temperature properties of QCD are on top of the list, since they mostly concern 
spatially homogeneous systems, where the equations of motion are still simple to solve. 
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1 INTRODUCTION 



There has been quite some understanding of the deconfmement transition in 4-dimensions 
on the basis of strings in strong coupling lattice QCD. It comes from a roughening of the 
strings due the entropic enhancement of configurations with long wiggly strings. The high 
temperature phase, however, is not understood in a picture where the underlying degrees 
of freedom are strings at low temperature. Indeed the old Hagedorn picture is limited 
to temperatures below the phase transition. Above the critical temperature T c spatio- 
temporal Wilson loops are no longer suppressed due to their area, the string tension goes 
to zero and therefore strings have seized to live in the plasma. There is a remnant of the 
low temperature theory in the behaviour of purely spatial Wilson loops, but this effect is 
not very strong at high temperatures. In Refs. [2] an analysis of the contribution due to 
these spatial surfaces has been made in 4- dimensions. There have been various attempts 
to analytically continue effective 4-dimensional string theories to the deconfmement phase, 
see e.g. Refs. [31 HJ E], however without any phenomenological application so far. 

The duality of string theory with 5-dimensional gravity can help. In conformal AdS^ 
the metric is well known. It has a horizon in the bulk space at = ^J- where (3 = 1/T, 
the inverse of temperature, and i is the radius of the AdS-space. Conformal solutions for 
entropy scale like s oc T 3 , since the 3-dimensional area of the horizon is given as A oc 2 r K 2 r\. 
Promising solutions of this conformal theory have been proposed to the problem of viscosity 
r] [6] with a small constant value for rj/s. Top to bottom approaches based on the conformal 
SYM with fermions have been investigating the chiral phase transition and problems at 
finite density [HE]. 

Indeed there may be a window with the plasma not close to the phase transition and 
not yet perturbative, where this conformal theory mimics truthfully reality. With solutions 
at T = at hand [H [9] which break conformal symmetry and give confinement, it is 
challenging to investigate the T ^ sector of the theory for all temperatures. Important 
progress has been made in this direction by two groups [THl [TT] . 

In this paper we follow this general approach, and investigate the equation of state with 
two different methods. On the one hand we use the Hawking-Page formalism [12J to derive 
the renormalized free energy from the action. On the other hand we directly calculate the 
entropy from the Bekenstein-Hawking formula [13J. These two approaches are assumed to 
give identical results as long as the black hole is treated entirely classical. Renormalization 
is well defined, because the T = theory has asymptotic freedom. We further investigate 
the change of the critical temperature with the number of flavours induced by the change 
of the running coupling constant. In a renormalization group framework interesting be- 
haviour of T c has been demonstrated recently [T4] ■ Finally using the string theory we can 
investigate the finite temperature behaviour of several thermodynamic quantities, like the 
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speed of sound, spatial string tension and vacuum expectation value of the Polyakov loop. 
Comparison with lattice data is made. 

The outline of the paper is as follows: In Section 2 we give an overview of the 5-d gravity 
action and specify the dilaton potential which is used. In Section [3] we will reproduce the 
zero temperature results focussing on the behaviour of the metric and the coupling near 
the boundary z — Y 0. Section H] gives the finite temperature calculation emphasizing the 
similarities and differences of the two solutions at the boundary. Section [5] is devoted to 
deduce the thermodynamics of the plasma using the Page-Hawking approach. In Section |6] 
we discuss the Bekenstein-Hawking approach to compute the free energy, and compare 
with the method of previous section. We also discuss the consequences of the flavour 
dependence of the running coupling. In Section [7] we compute the thermal behaviour of 
speed of sound, spatial string tension and vacuum expectation value of the Polyakov loop, 
and its comparison with available lattice data. Finally, Section[8]gives a final discussion and 
our conclusions. In Appendix A we discuss in details the procedure to solve numerically 
the Einstein equations at finite temperature. In Appendix B we show technical details for 
the analytical computation of thermodynamic quantities in the ultraviolet regime. 



2 5-d gravity action 

In the large N c limit, we assume a five dimensional gravity-dilaton model, with the action 

S = TtV f d " x ^ ( R - \ 9 ^ 9 ^ - V W)\ - f d 4 xVHK . (1) 

l07rG 5 J \ 3 J 87rG 5 J aM 

The last term is the Gibbons- Hawking term, where the integration is evaluated at 

the boundary dM of the five dimensional space given by z = 0. The induced metric on 

the surface is denoted by H. An infinitesimal distance z = e to the boundary is used 

to regularize our expressions. Later we will take the limit e — > 0. The metric G pu is 

taken in the Einstein frame, and the extrinsic curvature or second fundamental form of the 

boundary, K^ v , is evaluated with the help of the normal n p at the boundary dM: 



= -V^n v = \ p H^d a H^, (2) 
K = H^K, U . (3) 

In Ref. pQ we extrapolated the /3-function of QCD to the infrared with a parametriza- 
tion which was consistent with asymptotic freedom and the heavy qq potential at zero 
temperature. This parametrization takes the form 
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2 5-d gravity action 
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-b 2 a + 



b 2 a + ( — - 0o ) of 
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-a/a 



where 



a(z) 



(4) 



(5) 



is the running coupling. With this /3-function we could obtain the dilaton potential in 
AdS/QCD as a function of the running coupling constant: 



V(a) 



(3(a)\ \ fa\ — 



3a 



\a 



■Exp 
■Exp 



((2 7 -3)6 2 + 4/3 a + 2/3 1 a 2 ) 
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^e-S ( 36 2 - 4/3 « - 2/W + (- - 2f3 - 2^a)a 
9 \ a 



•Exp('^. El 



a 
ai 



(6) 



where E\ is the exponential integral function. For a < a the potential is strictly determined 
by the perturbative /3-function, 



13(a) = -f3 a 2 - f3 ia 3 - (3 2 a 4 + ... (7) 

The parameters /So and /3i are universal, i.e. they are regular izat ion scheme independent. 
For a > a the /3-function behaves linearly, [3(a) ~ — b 2 a, and the potential is characterized 
by the non-perturbative constants —b 2 and a. The parametrized coefficient (3 2 has the 
form 

6 2 - 3(3 a - 6/3i« 2 

P2 7—5 • l iU J 

ocr 

From the heavy qq potential at zero temperature one gets the optimum values [1] 

6 2 = 2.3, a = 0.45, 1 = 4.389 GeV" 1 . (11) 
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The potential itself approaches the conformal limit V = —12/£ 2 for a — > 0. Based on 
this action we will investigate the thermodynamics of QCD. 

Our parametrization given by Eq. (jlj) is simple and more intuitive than that of Ref. |15j . 
It allows for analytical computations in many cases, e.g. we can derive analytically the 
dilaton potential Eq. ([6]) from the /3-funtion Eq. (jl]). The dilaton potential of Ref. [15] is 
fine tuned and the role played by their parameters are not so obvious. We plot in Fig. [TJ 
the dilaton potential of Eq. (|6]) which we use in this work, and compare it with the one 
proposed in Ref. [15] which we call Vgkmn- We use in both cases the same value of I given 
by Eq. ( fill . Vgkmn(ch) is of the order of 10 4 GeV 2 in the regime of our interest, a ~ 0.3, 
in contrast to the value ~ 1 GeV 2 given by our potential V(a). Both potentials share the 
same ultraviolet behaviour, but they differ in the infrared. Since the scheme dependent 
coefficient 02 in Ref. [15] is very much larger than ours, the running of the coupling of 
Ref. [15] deviates already for extremely small values from the regime dictated by the leading 
coefficients Pq and fi\ and then it depends entirely on the "infrared" parametrisation in 
the regime of our interest. In the scheme of Ref. [15] the values of a sampling the dilaton 
potential V(a) are much smaller over the whole range of temperatures. At e.g. T = 5T C 
a(Ref. [15]) = 0.00095 or more than 100-times smaller than the usual running coupling in 
MS scheme a(7rT) = 0.11. This explains the large difference between the dilaton potentials 
when they are plotted as a function of the same a. Unfortunately no mapping between 
a(Ref. [15]) and a in our MS- scheme is known. A further discussion will be presented in 
Section El 



3 Thermal gas solution of the Einstein Equations 

The equations of motion corresponding to the five dimensional gravity-dilaton action, 
Eq. (Cp), are given by 



= T„ v , (12) 



with 



E/j,u — Rfj, u — -RG^ , (13) 

Tfiu = p^-^G^t^d^d^ + V^y (14) 

The left hand side of Eq. ( T121) is the Einstein tensor E^ v while the right hand side is 
the energy momentum tensor T^ u . The Thermal Gas solution preserves spatial rotational 
invariance and it has a metric similar to the zero temperature solution in Euclidean space: 
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3 Thermal gas solution of the Einstein Equations 



V(or) 




Fig. 1 : Dilaton potential V as a function of the running coupling a. The full (blue) line shows 
the dilaton potential V given by Eq. ^) with parameters Eqs. f3|), (TJ|) and 07]) . For comparison, 
we show as a dashed (red) line the dilaton potential Vgkmn given in Ref. fT5\/ multiplied by a 
factor 1CT 3 . 



G^u = ^(z)-diag (1,1, 1,1,1) , (15) 
ds 2 = bl(z) (dr 2 + dx k dx k + dz 2 ) , (16) 

where z G (0, oo) is the bulk coordinate in the fifth dimension, and the imaginary time 
coordinate r is periodic (r — > t + 0) with period = 1/T, the inverse of temperature. 
This solution exists at all temperatures. Under the assumption that the energy scale is 
proportional to bo, i.e. E = A^fro, the /3-function writes 

«") = «^. (17) 

There will be a second solution at finite temperature, the Black Hole solution, with a 
horizon in the bulk coordinate which characterizes the gluon plasma, and which we will 
discuss in the next section. After computing the Einstein tensor E^ v and the energy 
momentum tensor T^ u in terms of bo(z) and its derivatives, one ends up with three equations 
which determine the thermal gas solution: 



16 3 
W^ = -b W 2 + ~b V. 



7 



b' Q = -~b 2 W , (19) 
a' = a ^boWl. (20) 

In the following text ( ' ) stands for derivative with respect to the z coordinate. In 
order to reduce the equations of motion to first order equations, we have introduced the 
superpotential defined as 

W (z) = -l$\. (21) 

This definition agrees with the one given in Ref. pp. The functions Wo, b and 
characterizing this solution have the index 0. We replace the variable <ft, dilaton potential, 
by the running coupling a . 

The system Eqs. (|18|) - (|20l) has been solved numerically in Ref. [T] by considering the 
bulk coordinate at = 0.3426 GeV -1 which is mapped to the energy E* = 3 GeV by the 
metric factor 6 an d the arbitrary scale = 0.264 GeV. The choice of scale size A# has 
only historical reasons, and in principle the mapping of the energy coordinate to the z- 
coordinate has an arbitrary constant. The resulting parametrization of the /3-function and 
dilaton potential would then also change. We use this scale, because we did not want to 
recalculate the fit to the running coupling and the string tension with another energy scale. 
The value of the running coupling at E* = 3 GeV follows from the experimental PDG data 
of the running coupling, c.f. Ref. [16], and it is ao(z*) = 0.25241. The value of Wq comes 
from the Einstein equations which allow an explicit solution of the superpotential in terms 
of V(a) and /3(a), 



b (z*) = E*/A E , (22) 
a (z*) = 0.25241, (23) 

We have solved the above equations and verified that they give the same solutions as 
obtained in the previous paper on zero temperature Ref. [1J. Note the finite temperature 
— 1/T does not enter in the gravity equations dependent on z. The thermal gas is solely 
defined by the periodicity in r. The initial conditions fix the running of the coupling in the 
bulk and its scale A which we compute from the numerical solution. Given the perturbative 
/3-function, which is parametrized for higher orders in Ref. [T] 

(3(a) = -(3 a 2 - fca 3 - f3 2 a 4 + ... , (25) 
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4 Black hole solution of the Einstein Equations 



the running coupling has the form: 



a Q (z) 



k + (-i log(L H^ + 



+ (|(log(L z )) 2 - + 2/3^) log(L 2 ) + c^j ± + 0(L: 



(26) 



with 



^ +ftfc2+ i*-i + !' (27) 

and the definition 

L 2 :=-log(zA). (28) 

The ultraviolet expansion of bo(z) follows from Eqs. (|T9|) and (1201) . using for «o the 
expansion given by Eq. (1261) . Then 6 reads 



-/W^) + (22/3 2 - 9ft) a 2 ^) 



(602/3 3 - 540/3 /3i + 81/3 2 ) a 3 (z) + 0{ 



2187 



On 



(29) 



In order to compute accurately the value of A, one can consider b as a function of 
ao(^) and expand it around the point cto(-2*) = a* where the initial conditions were given. 
The expansion writes 



6 («o) 



Est 



1 f . (1 -£'(«.)) 

■(a Q - a* J + 



+ 



/3(a*)"~ u ' 2(/3'(a,)) 2 

1 - 3/3'(a*) + 2(/3'("*) - /3(a*)/3>*)) 2 ) 



(a - a*) 2 
(at - «*) 3 + ■ 



(30) 



6(/3(a*)) 2 

For Nf = a fit of the numerical solutions for aco(z) and bo(z) to the form given in Eqs. ( |26l) 
and (|30|) respectively, yields in both cases A = 0.543 GeV which is a factor two larger than 
l/£ — 0.228 GeV. The parameter k is an integration constant which appears in the Gell- 
Mann-Low integral [17] . Note that a particular choice of k fixes the definition of A. We 
have chosen in this computation k = 0. 



4 Black hole solution of the Einstein Equations 



The gravitational equations have two different types of solutions for the metric. Besides 
the thermal gas solution which we discussed in the previous section, there is also a solution 
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which has a horizon localized in the bulk coordinate at z = Zh similar to the situation in 
4-dim gravity. The phenomenology of the gluon plasma arises from the competition of the 
free energies computed in both metrics. When the free energy of the black hole solution 
has a lower value than the thermal gas solution, the phase transition to the quark gluon 
plasma takes place [1~2~] . 

Now we will discuss the equations of motion at finite temperature using the black hole 
metric. The procedure is similar as discussed above for the zero temperature case, but a 
new equation appears due to the black hole factor f(z). The black hole metric in Einstein 
frame has the form: 

Gfx V = & 2 (z)-diag(/(z),l,l,l,^j , (31) 

ds 2 = b 2 (z)(f(z)dr 2 + dx k dx k + ^-], (32) 
V /(«)/ 



where 



/(0) =1, (33) 
f(z h ) =0. (34) 



Near the horizon the metric is given by 

dz 



ds 2 = b 2 (z h ) f'(z h ) ■ (z - z h )dr 2 + dx ■ dx + — — ) . (35) 

V / \ z h) ■ [z - Zh),' 

We define a new variable p := y/z — Zh- In terms of p we obtain 
2 46 2 (^) / 2 (f{z h ), x , ,, ... 



ds 1 = - ; 7 p 2 J -^^dr + dp 2 + b 2 (z h )dx ■ dx. (36) 
f'(zh) \ V 2 J 

The r-p portion of the metric defines a two-plane in polar coordinates with r serving 
as the angular coordinate. To avoid a conical singularity at p = we must require that 
\f'(zh)T /2\ has a period of 2n. In Matsubara (imaginary-time) formalism the period is 
equal to inverse temperature /3 = T -1 . Thus, the temperature of a black hole solution is 
given by PHCL9] 

r = _Z^l. (37) 

The Einstein tensor E^ u and energy momentum tensor T^ u have different components 
in 00, 44 and spatial directions, and can be expressed in terms of b(z), f(z), dilaton field 
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4 Black hole solution of the Einstein Equations 



<f>(z), dilaton potential V((p) and its derivatives. One ends up with the four equations 
relevant at T ^ 



W = ^bW 2 - j (wf - \bV^ , (38) 

b' = - -b 2 W, (39) 
9 

a' = aVbW', (40) 

f" = p'bW , (41) 

where we have introduced the superpotential at finite temperature which is defined in 
analogy with Eq. ( )2T|) . 

(42) 

Note that the system of Eqs. (!38 l) -( l4B reduces to the zero temperature formulas 
Eqs. (|T8l)-(l20l when / = 1. Our prescription is to use the same dilaton potential V(<f>) 
at zero and finite temperature. Note, however, that V(z) is affected by the tempera- 
ture dependence of a(z) = e*W, i.e. V(z) = V(a(z)). In order to solve the system of 
equations (!38|) -( l4T]) one should specify five initial conditions, as one has to handle three 
differential equations of first order, Eqs. ()38l)-(l40l. and one differential equation Eq. ( )4T1) 
of second order. The Eq. (jUJ) for the black hole function f(z) can be solved analytically 
in terms of b. Since our set up includes asymptotic freedom, b(z) is different from the 
conformal solution b coni (z) = Ijz. Also the black hole factor deviates from the simple form 
of the conformal solution f cont (z) = 1 — z 4 / ' z\. Using Eq. ( 139]) one gets 

J o y u > Jo 

where the first integration constant Cf has been chosen such that f(zh) = 0. The second 
integration constant is fixed to unity by the requirement /(0) = 1. From Eq. ( )37|) and the 
definition of Cf given by Eq. ( 143|) . one derives easily the useful relation 

C f = AixTb\z h ). (44) 

The technical procedure to solve numerically the Einstein equations is explained in 
details in Appendix A. We show in Figs. 121 El S]and|5]the solutions for W(z), z ■ b(z)/£, 
a(z) and f(z) obtained at the temperature T = 368 MeV and compare them in the same 
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Fig. 2: Superpotential W as a function of z. The full (blue) line gives the finite temperature result 
at T = 368 Me V, and the dashed (red) line the zero temperature superpotential. The maximum 
value of z shown in the plot corresponds to the horizon Zh(T) = 0.946 GeV^ 1 . 

figures with the zero temperature solutions as computed in Sec. El Note that the T = 
and T ^ solutions agree in the ultraviolet, i.e. z < 0.5 GeV -1 , but differ in the infrared 
due to the thermal fluctuations introduced by the black hole horizon at Zh = 0.946 GeV -1 . 

As we will explain in Sec. El we use a(z) to extract the trace anomaly. The finite 
temperature solution differs from the zero temperature one at order 0(z A ). Care has to 
be taken to keep track of leading logarithmic effects which are not usually considered, 
c.f. Ref. [10J, but are important if one wants to calculate the gluon condensate G. The 
difference between zero and finite temperature solutions is mainly given by the gluon 
condensate or the trace anomaly which equals G up to normalization factors, 



G 



7t(j5 (3(a) 



15 



a" 



{(TrFp T - (TrF^)o) . 



(45) 



The subscripts T and stand for the thermal and vacuum values of (TrF^ u ) respectively. 
The expressions relating zero and finite temperature quantities read up to 0(z 4 al) 



b(z) 
a(z) 
W(z) 



b (z) 
ato(z) 
W (z 



l + ^z 4 (l + c b 1 a (z) + c b 2 a 2 o(z)) + 
45 G 



'P/3 a (z) 



a (z) + c 2 a {z) + 



1 



5(7 1 f l + c\a Q {z) + U + ^ - l ^Bl ) al(z) ) I 



45' 



(46) 
(47) 
(48) 
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4 Black hole solution of the Einstein Equations 




Fig. 4: Running coupling a as a function of z. We show the finite temperature result at T = 
368 Me V and the zero temperature one. See Fig. [D for convention. 
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Fig. 5: f as a function of z. We show the finite temperature result at T = 368 MeV. 



z^QAz) 1 - If z " f 1 + c i«° W + ( c 2 - T A + ?) <*)) 
f( z ) = l ~ aA (A ~~r^ 7 ; h + ( 49 ) 



where 



9 „ ft\ 4 2 ft ft 2 ft 



being Qf(z) given by (c.f. Appendix B) 

Qf(z) = 1 + ^/3oaoW - ^(7/3 2 - 6ft)ag(z) + O(a 3 ) . (51) 

Higher orders in «o and 2 in Eqs. fH6l)-(l49l) are indicated by dots. 

To arrive at Eqs. ( |47|) -( |49|) we have assumed for b(z) an expansion of the form given 
by Eq. ( H6|) . and used the equations of motion Eqs. ( l39]) - fpf0l) . Note that on the r.h.s. of 
these expressions, inside the brackets, one can substitute a (z) by a(z) and the expressions 
remain valid at this order, as the difference between both quantities, cto(z) and a(z), is 
0(z 4 ), c.f. Eq. (H7|) . To get the values of the coefficients c\ and c\, one has to substitute the 
expansions for b(z), W(z) and f(z) into the first equation of motion, Eq. fl38|) . and use the 
assumption that the dilaton potential at finite temperature as a function of the dilaton field 
has the same functional form as the one at zero temperature, i.e. Vr(a) = Vr=o(«o)|a =a- 
Then one gets 

4 = j^A. (52) 
4 = -^ + I ft -jSL/e, (53) 
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5 Free energy from the Einstein-Hilbert action 



where Cf is defined in Eq. (l4"3j) . Note that the expression of c\ in Eq. (|53|) means that in 
the expansion of the quantities in powers of z, Eqs. P6|) - (|49|) . there are contributions not 
only of the gluon condensate G, but also of Cf. The leading contribution involving Cf is 
~ z 4 Cfa.Q(z) times the corresponding quantity at zero temperature (~ z 4 Cfao(z) in the 
case of a(z)). 

The contribution of G in Eqs. f|46|) - fj49|) is visible in Figs. [2]|H and all these figures 
consistently show that G is positive. Note that in our set-up the correction arising from the 
NLO and NNLO coeficients of the z 4 -term are not small in the infrared, since oeo(zh) ~ 0.5 
at T = 368 MeV. Therefore we have to resort to a UV-analysis to determine the gluon 
condensate G from the computation of a(z). 



5 Free energy from the Einstein-Hilbert action 

To get the thermodynamics of the five dimensional gravity-dilaton model of Eq. ([Q) one 
computes the free energy at fixed temperature by introducing a lower cut off z = e in the 
integral over the on shell action: 

&T = S ICS (e) . (54) 

Regularization is needed due to ultraviolet divergences near the holographic boundary 
e — > 0. The procedure to compute the regularized action with the black hole solution is 
explained in details in Ref . [10] ■ The free energy is computed as the difference between the 
free energy of the black hole solution and that of the thermal gas solution, so by definition 
the later has zero free energy. The result for the free energy is [10] 

T = \ BmraW - «W) - (l6G - °i) • (55) 

From J 7 one can compute the pressure and the rest of thermodynamic quantities by 
applying the thermodynamic relations. The values of Cf and G in the ultraviolet are given 
by (cf. Appendix B) 



G 



4£ 3 



4vr¥ ! T 4 



p a h + i(ll/3 2 - + 0(a 3 h ) 



71 



l £ 3 



45 



rp4 



(5la 2 h + 0{ai) 



(56) 
(57) 



Inserting Eqs. (1561) and flBTl) into Eq. f[53|) one gets the UV assymptotic expansion of the 
free energy which will be discussed later. To deal with Eq. ( 155]) at temperatures near T c 
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one has to compute the temperature dependence of Cf and G numerically. Using Eq. (jHJ) 
one can calculate Cf from the numerical result of b(z). On the other hand the temperature 
dependence of G can be computed by comparing the thermal gas and black hole solutions 
in the ultraviolet, using Eqs. (I46l) -( j49l) . In the UV we perform a fit of the difference 
a(z) — Qq(z), and compute the coefficient G for different temperatures. Note that it is 
much more efficient to use a instead of b, as the latter diverges in the UV making it more 
difficult to get reliable results for G. In this paper we have analyzed carefully the expansion 
of the 0(z 4 ) term in a(z) —a>o(z), Eq. ( 1471) . which gives the trace anomaly G in the plasma. 
Higher order terms in «o affect appreciably the fit of G, and it is indispensable to consider 
at least the order 0{olq) to get good agreement of the thermodynamic quantities with the 
numerical results from the Bekenstein-Hawking entropy formula as computed in Sec. [6j Up 
to C(«o) it reads 

a(z)-a (z) 45 G ( ( b B (3A , . 2f A 

— ? — " Tm v v t~ aJ ao{z) C2a ° {z) ) ' ( } 

where c\ and are given by Eqs. (|52p and (15171) respectively. In Fig. [H] we show the 
numerical and analytical results of (oc(z) — ao(z))/z A for small z. Note that this quantity 
is not flat in this region, and a rough fit with a constant term ~ G (constant in z) leads 
in general to an overestimation of the value of G, and then also on the value of J 7 , cf. 
Eq. ( )55l) . This induces an appreciable error in the behaviour of J 7 , and in the value of T c . 
For instance, if we performed the fit by neglecting all higher orders in cto on the r.h.s. 
of Eq. fl58|) . we would get for the transition temperature T c = 298.7 MeV, instead of the 
correct value T c = 273.0 MeV for zero flavours. 



6 Free Energy from the Bekenstein-Hawking Entropy 

One of the postulates of the gauge/string duality is that the entropy of the hot gauge 
theory equals the Bekenstein-Hawking entropy of its gravity dual. This opens a quick way 
to check the results of the previous section independently. The Bekenstein-Hawking entropy 
is proportional to the (3-dimensional) area of the black hole at the horizon r = = £ 2 / Zh 
with the metric Eq. ( 13TT) . In non conformal AdS/QCD the entropy writes 

S(T) = ^& 3 (r,(T)) , (59) 
where G5 is the gravitational constant in 5-dim and b is the metric factor in Einstein frame. 

2{' 

The pref actor \/{4G^) = in conformal theory is much too large for pure QCD, since 

Af = 4 supergravity includes extra degrees of freedom, gluinos and scalars, besides gluons. 
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Fig. 6: Difference between the running coupling at finite and zero temperature divided by z A . 
We consider T = 368 Me V. Blue line corresponds to the numerical computation from Sees. [3] 
and R[ We plot as a dashed blue line the regime z < 0.02 GeV~ 1 which is affected by numerical 
errors. Dashed red line corresponds to the r.h.s. of Eq. I158\) including up to 0(oq). We also 
show as a dashed dotted black line the r.h.s. of Eq. A58\) . but neglecting all the orders in ao, i.e. 
-45G/(8£ 3 /3 ). 



In order to map out an equation of state, one needs the location of the horizon = £ 2 / Zh 
as a function of temperature. One finds for a temperature above some minimum value T min 
in general a solution with a small (small black hole) and a solution with a large 
horizon (large black hole). Only the large black hole is stable, because its free energy is a 
minimum. The large black hole solutions are used to calculate the entropy. 

The free energy due to black holes must be calculated by combining the entropy Sb of 
big black holes and the entropy S$ of small black holes. The free energy of the big black 
hole can then be computed as [10J 



T B = Moo)- S s dT- S B dT, (60) 

where the unstable free energy from small black holes vanishes in the limit T — > oo, i.e. 
Js(oc) = 0. 

The calculation of the free energy in terms of the entropy is much easier than the full 
calculation shown in section [5j The reason is that the extraction of the gluon condensate 
G is rather subtle and we had to develop the procedure indicated in the previous section to 
get a reliable value for G. The perfect agreement of both methods gives a good guarantee 
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Fig. 7: Free energy density as a function of temperature. We show as a full (red) line the 
result obtained starting from the Bekenstein- Hawking entropy formula, c.f. Sec. Blue points 
correspond to the result using the Einstein- Hilbert action, c.f. Sec. O Eq. Ii55\) . We include up 
to 0((Xq) in the r.h.s. of Eq. [5B\) to compute G. 

that the numerical solutions are reliable. 

We show in Fig. [7] the free energy obtained by using the Bekenstein-Hawking entropy 
formula (full red line) as computed in this section, and the Einstein-Hilbert action (blue 
points) from Section [5j using NNLO terms in a in Eq. ( |58l) to compute the gluon conden- 
sate G as a function of temperature. One clearly recognizes in this figure the first order 
phase transition at the temperature T c = 273 MeV for zero flavours, which is quite close 
to lattice simulations. We consider the scale setting from the zero temperature gravity as 
a great success for gauge gravity duality. The upper branch in Fig. [7] represents the small 
black holes which are energetically disfavoured. 

For high temperatures a weak coupling expansion of the pressure, p = — J 7 / Vol (3), can 
be made. We refer the reader to Appendix B for a detailed discussion of the ultraviolet 
properties of the thermodynamic quantities. For an analytical computation in the ultravi- 
olet, one makes a weak coupling expansion in at = cx(zh)- The coupling at evaluated at 
the black hole horizon is finite. One can relate ath with the value of the running coupling 
at the scale z = l/{nT) by the following equation: 



where «t is defined as indicated. The easiest way to compute the pressure p(T) from the 




(61) 
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entropy density s(T) is to solve the equation dp(T)/dT = s(T). Then one can consider a 
general scheme for the weak coupling expansion in the pressure, p(T)/T 4 = p + picth + 
P2&h + + ' ' ' > an d solve for the coefficients in the above equation with s(T) given by 
Eq. (11091) . dp(T)/T can be computed easily by making use of Eq. (II 111) . Then one can 
identify the coefficients in the expansion. The result is 



p(T) 7T 3 £ 3 
~T r 16G 5 



1 - poot h + I (4/3 2 - 3ft) a\ - ^ (91/3 3 - lUfoA + 72/3 2 ) a\ + O(o*) 
[l - 2.33a h + \Ma 2 h + 0.33aj[ + 0(a£)] , (62) 



where we show in the last expression the values of the coefficients corresponding to N c = 3 
and Nf = 0. 

In AdS/QCD we may choose the gravitational constant to reproduce the ideal gas limit 
at high temperatures, i.e. p(T)/T 4 ~t->oo — l)7r 2 /45. Then one gets 

(63) 



16Gg° 45tt£ 3 

We call the so determined constant Gf^ indicating that it follows from the large temper- 
ature limit. As pointed out in Ref. this value is a factor 8/45 smaller than the value 
for M = 4 super Yang-Mills theory. This decrease may be explained by the following two 
arguments. The number of degrees of freedom is reduced in QCD compared with SQCD by 
a factor (2/15). QCD is weakly interacting at high energies compared with the AdS/CFT 
theory which remains strongly interacting. This gives another factor (4/3). 

We show in Fig. [8] the pressure as a function of temperature, normalized to the Stefan- 
Boltzmann limit. It is noteworthy and visible in the figures that the expansion in terms 
of ah converge quite rapidly. On the other hand, as one can see the holographic model 
approaches the ultraviolet limit slower than lattice data. The asymptotic expansion of the 
pressure in QCD with iV c = 3, Nf = has the form: 



Pqcd{T) _ 8tt 2 



15 on /a\3/2 
— a + 30 - +... 

47T V 7T / 



, . . . . (64) 

If one compares the 0(a) coefficient in the ultraviolet expansion of the holographic 
model (pi dS = ||) Eq. (1621) with the corresponding coefficient from perturbative QCD 
(pf co = g) ? one gets a factor pf s /pl QCD = 88/45 ~ 1.956. This ratio explains the 
deviation observed in Fig. [HI since the leading coefficient gives a good approximation to 
lattice QCD at high temperatures T = (100 — 1000)T C . The same factor appears when one 
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Fig. 8: Pressure normalized to the Stefan- Boltzmann limit, as a function ofT/T c . We show as 
points the high temperature lattice data for SU(3) taken from Ref. [20]. We plot the analytical 
result from the holographic model for several orders as an ultraviolet expansion in powers of the 
running coupling, c.f. Eq. h62\). and also the QCD perturbative result up to order 0(a), c.f. 
Eq. (6%). 



compares the first non-zero coefficient in the ultraviolet expansion of the trace anomaly 
which is 0(a 2 ), and so one expects that the holographic model also predicts values of 
(e — 3P)/T 4 larger than lattice data at high temperatures. This seems to be a property of 
the present holographic model, and there is no easy way to cure it. In Ref. [TT] a behaviour 
/3(a) ~ —a q for q > 2 has been studied and numerical consistency between QCD and this 
AdS-model with a very simplified /3 function is reached for q = 10/3. This parametrization 
disagrees, however, with the standard perturbative behavior of the /3-function of QCD, 
which has been the basic starting point to make AdS similar to QCD in the program of 
Refs. [SHEaCE]. 

The parameterization presented in Ref. [15] is successful to reproduce lattice data in 
the regime 1 < T/T c < 5, but in view of our analysis it is clear that this is because non 
perturbative effects are much stronger than perturbative ones even at very high tempera- 
tures, which seems to be not reasonable. As a matter of fact, the /3-function of Ref. [15] 
agrees only for a s < 10~ 7 with an 1% accuracy with the asymptotic /3-function of QCD. 
The dilaton potential used in this reference has the correct uv-behaviour of QCD, but for 
all practical purposes it is a model potential which is designed to fit the thermal equation 
of the gluon plasma and not the /3-function established in perturbative QCD for a s < 0.2. 
Since in Ref. [15] the scheme dependent /?2« 4 term in the /3-function is very large, it is 
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simply not clear how to determine the running of a s in the MOM or MS scheme from the 
parametrization of /3(a) in this reference. 

It is interesting to analyze the dependence of the phase transition temperature when 
the number of flavors is changed. The authors of Ref. [H] have studied the iVy-dependence 
of the transition temperature T c with the help of renormalization group flow equations. As 
shown in this reference, the scale A QCD changes with Nf. Therefore it is recommended to 
solve the Einstein equations by keeping the running coupling fixed at a mid scale 3 GeV, 
c.f. Eq. ( 123]) . In order to study the flavour dependence of T c in the present model, we vary 
the number of flavors Nf from Nf = to Nf = 10 in the coefficients (3 and (3i of the /3- 
function, c.f. Eqs. which control the short distance (high temperature) regime. We 

retain the nonperturbative parameters 62 and a as in Eq. (1111) . which is reasonable because 
62 and a are responsible for the infrared large distance (low temperatue) regime, and the 
string tension in the qq potential mainly depends on these two parameters. We don't study 
the effect of dynamical quarks, therefore our analysis is restricted to a quenched approach. 

Fig. [9] shows the dependence of the critical temperature T c on the number of flavors 
Nf < 10. Larger values of Nf are difficult to implement numerically. We get for Nf = 
as transition temperature 

T Nf=0 = 273.0 MeV, (65) 

which is very close to the lattice results T c = 270(2) MeV [23]. It is gratifying that the 
absolute value of the transition temperature comes out so well inspite of the slow con- 
vergence of the pressure towards the Stefan-Boltzmann limit with increasing temperature. 
Comparing results with different flavour numbers we obtain an almost linear behaviour of 
the transition temperature for small Nf-. 

T c = T Nf=0 (l - nNf + 0(N 2 f )) , k = 0.1205 . (66) 

This linear scaling of the critical temperature with Nf for small Nf has been claimed in 
Refs. [TH [24]. The value of k we get is quite close to the one estimated in Ref. [14] , 
K ~ 0.107. The flattening of the function T c (Nf) for high values of Nf in Fig. |9] is in 
accordance with Ref. [H]. This reference explains this fact as a consequence of the IR 
fixed-point structure of the theory. In our case the /3-function given by Eq. (j3J) does not 
have an IR fixed-point, and the flattening is a consequence of the weakening of the infrared 
coupling a when Nf increases. This means that it takes smaller scales to reach a critical 
coupling a(T c ) to bind the gluons into glueballs, and therefore T c has to decrease. We 
found that for the points plotted in Fig. [9j the value of a(T c ) is not affected much by Nf. 
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Fig. 9: Phase transition temperature as a function of the number of flavors. The points correspond 
to the numerical computation of the holographic model as explained in Sees. [3| and El within a 
quenched approximation. We also plot as (red) continuous line the scaling law of Eq. 



7 Thermodynamic Observables independent of G : 



We study in this section several thermodynamic quantities which are independent of the 
5-d gravitational constant G§. For a more complete discussion of those thermodynamic 
quantities which dependent on G5 we refer to Ref. [25] . 

First we focus on the speed of sound c s . From the specific heat per unit volume c v : 



T 



d 2 p 
df 



2 ■ 



(67) 



and the entropy density s, one obtains the speed of sound: 



s 



(6S 



A computation of this quantity in the ultraviolet leads to (see Appendix B for details) 

(69) 



s 



\pWh + l(3o {Pi - 4/?i) al + 0(ai) 



In Fig. HO] we show the speed of sound computed with the holographic model, and 
compare the result with lattice data of Ref. [26] . We also plot the analytical ultraviolet 
approximation given in Eq. (1691) including several orders in an expansion in oth- Since (? s 
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Fig. 10: Speed of sound squared as a function of T jT c . We show as points the lattice data for 
SU(3) taken from Ref. [26]. The colored curves represent the analytical result from the holographic 
model for several orders as an ultraviolet expansion in powers of the running coupling, and the 
black solid line refers to the full numerical result. 



becomes close to 1/3 in the calculation, we see that we have massless excitations in the 
plasma in the temperature range 2T C < T < 5T C . 

The spatial string tension is another quantity which is very useful to test AdS/QCD 
models. It is non-vanishing even in the deconfined phase, and it gives useful information 
about the non perturbative features of high temperature QCD. With a quark and an 
antiquark located at x = | and x — — | respectively, the computation of the correlation 
function of rectangular Wilson loops in the (x, y) plane leads to a potential between quark 
and antiquarks which behaves linearly at large distances, i.e. 

(W[C\) y ~°° e-y y ^ , V{d) a s -d. (70) 

For details on the computation we refer the reader to e.g. Ref. [27] and references therein. 
The spatial string tension takes the following form: 

<?) = ^Tn*H) , (ri) 

where l s is the string length. Making use of all the technology developed in Appendix B, 
we can easily compute the UV asymptotics of o~ s (T). Using the ultraviolet expansion of 
b(a) given by Eq. f)102p and the corresponding expansion of given by Eq. fllOSp . then 
Eq. dZTjl leads to 
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1:2 



4 



2187 



(931$ - 1836AA + 642/3 2 ) a\ + 0{a\) 



(72) 



We show in Fig. [TT] a plot of T j ^/d~ s as a function of temperature including several orders 
in Eq. f!72|) . and the full numerical computation from Eq. f JTTj) . We show for comparison 
also the numerical result obtained from the model of Ref. [27]. One can see that our 
model reproduces very well the lattice data in the regime 1.10 < T/T c < 4.5. A fit to the 
lattice data taken from Ref. [22] gives a good xVd.o.f. < 1, and it is obtained by using 
l s = 1.94 GeV -1 which is 30% larger than the value quoted in Ref. [1] based on a joint 
analysis of the heavy QQ potential and running coupling at zero temperature. From a 
computation of the string tension at zero temperature one can see that this increase in l s 
can be partially explained as an effect of the change in the number of flavors, as Ref. [1] 
considers Nf = 4 while we consider here Nf = 0. The string tension at T = can be 
computed as [27] 

o = -^bl{z*)ocT{ z *) , (73) 

where gives the minimum of &o(,2)a: 4//3 (z). Using the solution of bo and otQ for Nf = 0, 
one reproduces the physical value o ~ (0.420 GeV) 2 for l^ f ~° = 2.22 GeV -1 , while for 
Nf = 4 one gets lf f ~ 4 = 1.45 GeV -1 pQ. The discrepancy of l s with the value we get 
from the fit of a s (T) is 12%. 

The vacuum expectation value of the Polyakov loop serves as an order parameter for 
the deconfinement transition in gluodynamics. The correlation function of two Polyakov 
loops taken in the large distance limit leads to the vacuum expectation value of one single 
Polyakov loop squared. This means that the Polyakov loop is related to the free energy of 
a single quark F q as 

(V(x)} = e~^ ] . (74) 

One of the main problems of the computation of this quantity is the renormalization. 
The multiplicative renormalizability of the Polyakov loop was first established in Ref. [28] . 
The Polyakov loop was computed in perturbation theory for the first time in Ref. [29] 
in pure gluodynamics. There has been recent progress to renormalize it in the lattice 
following different methods based on the computation of the one point and two point 
correlation functions of Polyakov loops. The multiplicative renormalization is then reached 
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Fig. 1 1 : T '/ ' ^fo~ s as a function of temperature (in units of T c ). We show as filled (red) points 
the lattice data for SU(3) taken from Ref. [26]. The colored curves represent the analytical result 
from the holographic model for several orders as an ultraviolet expansion in powers of the running 
coupling, c.f. Eq. \12^ , and the black solid line refers to the full numerical result of Eq. |7ip . We 
use the value l s = 1.94 GeV^ 1 . We show for comparison as square points the numerical result 
obtained from the model of Ref. JF7[ /. 



by identifying and extracting the quark self energy, see e.g. Refs. [301 EH]- in Ref- [22] it 
was proposed by one of the authors a phenomenological ansatz based on a dimension two 
gluon condensate of the dimensionally reduced effective theory of QCD, which was quite 
successful to reproduce lattice data of the Polyakov loop in the deconfined phase down to 
T = 1.03T C . This ansatz follows from the observation for the first time in Ref. [32] that 
close and above T c the behavior of the Polyakov loop is characterized by power corrections 
in 1/T 2 . These power corrections were later observed also in the equation of state of 
gluodynamics [35]. The computation of the Polyakov loop within the AdS/QCD formalism 
was addressed recently in Ref. [36] within a model based on a specific choice of the warp 
factor b(z) which naturally introduces these power corrections. In the following we will 
consider this approach, but using our model dictated by the 5-d grativy action. 

One can compute the vacuum expectation value of the Polyakov loop from the Nambu- 
Goto action of a string hanging down from a static quark on the boundary into the bulk. § 

1 See also Refs. |33II34) for an application of this ansatz to compute the heavy quark- ant iquark free energy 
and the equation of state of QCD. 

2 The coupling of the two dimensional curvature R 2 to the dilaton field is a well known a' correction to 
the Polyakov action which propagates to the Nambu Goto action [37]. In first order it enters there by a 
modified non conformal metric, as it is considered here. Higher order terms are interesting to investigate, 
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The fundamental string is stretched between the test quark at the boundary [z = 0) and 
the horizon (z = z^) of the black hole solution. See Ref. [36] for details. The Nambu-Goto 
action then reads 

1 r Zh 

S " G= 2^Tj ^« 4/ W(2)- (75) 

The action Eq. (j75p is divergent at z = 0. One can regularize it by substracting the 
action of the thermal gas solution up to a cutoff z c : 



Orcg 
NG 



1 



2nl 2 s T 
1 



dz a 4/ \z)b 2 {z) - / dz ai /3 (z)bl(z) 



LJO 



dz (a i/3 (z)b 2 (z) - a /3 (z)b 2 (z) 



LJO 



dz a^ S \z)1)q{z) 



(76) 



2nl 2 s T 

The cutoff becomes necessary because the free energy of a single quark diverges at 



T = 0. Note that z c introduces a normalization constant into the free energy F q = T ■ S : 



NG' 



In the second equality of Eq. (!76|) we have divided the action corresponding to the thermal 
gas solution into two integrals. The first integral inside the bracket in Eq. (1761) is UV 
convergent, as zero and finite temperature solutions have the same behavior in the UV. 
The renormalized vacuum expectation value of the Polyakov loop then writes 



L R (T) = e-^o . (77) 
S^q defined in Eq. (1T6"|) tends to zero in the limit T — )• oo independently of the value of 



z c , and so Lr(T) tends to 1. We show in Figure [15] as a continuous black line the behavior 
of Lr as a function of temperature computed numerically from Eqs. (1761) - (1771) . and its 
comparison with lattice data for gluodynamics with N c = 3 taken from Ref. [3l|. In order 
to reproduce lattice data, we have performed a fit by considering the string length l s and 
the cutoff z c as free parameters. The best fit in the regime T c < T < 10T C leads to 



L = 2.36 GeV~ 



0.43 GeV~ 



(78) 



Note that l s is 20% larger than the value we used for the spatial string tension, but it differs 
just 5% from the value one needs to reproduce the string tension at T = (see Eq. (!73|) 
and discussion below). Our approach fits the Polyakov loop very well without a dimension 
two condensate, since a dimension two operator would show up in the ultraviolet expansion 
of the thermal solutions near z = 0, Eqs. (H6l) - (H9l) . This does not exclude that a good 



but are in the scope of a separate longer study. 
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Fig. 12: Expectation value of the Polyakov loop as a function of temperature (in units ofT c ). 
Full (black) line corresponds to the numerical computation of Eqs. J 76] }- (D). We show as points 
lattice data for SU(3) taken from Ref. J22F for N% x N T = 32 2 x 4 and 32 3 x 8. The colored 
curves represent the analytical result from the holographic model including several orders in a^, 
c.f. Eq. fl79[ ). We also plot the result from standard perturbative QCD up to 0(a^) given by 
Eq. fEQ). We use the values l s = 2.36 GeV' 1 and z c = 0.43 GeV 1 . 



fit to the data exists of the form — 21ogL# ~ a + b(T c /T) 2 with a = —0.23, b = 1.60, in 
accordance with Ref. [32] (see also Ref. [36]). The Polyakov loop is zero in the confined 
phase and our approach gives a nonzero value at T c given by Lr(T c ) = e~5>( a+fe ) = 0.50. 
This first order jump is similar to the one predicted by the more reliable lattice data 
x N T = 32 3 x 8, c.f. Fig. M 

We compute in details in Appendix B the UV asymptotics of the Polyakov loop. The 
results is 



L R (T) = exp 
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a? 

2 h 



1 + ^Poa h + ^ (l61/3 2 + 72ft) a\ + 0(a 3 h ] 



(79) 



We show in Fig. [12] the analytical result given by Eq. (JT9"j) . using the value of l s quoted 
in Eq. (|78l) . The Polyakov loop was computed in perturbative QCD up to NLO in Ref. [29] 
(see also Ref. [22]), and it has been recently corrected by two groups Ref. [Ml 39J. For 
gluodynamics this gives 



L PT (T) = exp 



2N r . 



3 

n <\f + 



Nl-l 



log «t + log(47r) 
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+ 0(a 



5/2 x 
T ) 



10) 



27 



Note that since the perturbative /3-function starts at order a 2 , changes in the scale 
/i affect 0(a 5 / 2 ). In Eq. ([79]) the power counting in doesn't follow the perturbative 
scheme. This discrepancy with PT is common of all the renormalization group revised 
models constructed by the general procedure of Kiritsis et al., c.f. Refs. [22] EI]- We 
have plotted in Fig. [12] also the perturbative result given by Eq. flHUj) . Note that lattice 
data approach the perturbative result very accurately for T above 8T C . Here the AdS- 
perturbation theory does not seem to converge rather rapidly. 

8 Discussion and final remark 

We have demonstrated in this work the numerical agreement between computations of the 
free energy from the Einstein-Hilbert action and from the Bekenstein-Hawking entropy 
formula. Both approaches leads to the same result, but the former method is much more 
sensitive to numerical errors, and an accurate computation is only possible when one takes 
care of including leading logarithmic effects in an ultraviolet expansion of the scale factor 
at finite temperature. 

We have also computed analytical expressions in the ultraviolet for the thermodynamic 
quantities as an expansion in powers of the running coupling ah evaluated at the black 
hole horizon. This expansion turns out to converge quite rapidly even at temperatures 
T ~ 1.5T C , quite opposite to the conventional QCD perturbation theory at high tempera- 
ture [ID]. We have extended our analysis to other thermodynamic quantities computed in 
the string frame, in particular the spatial string tension and the vacuum expectation value 
of the Polyakov loop, and the agreement with lattice data is better in this case. 

From our analysis we see that the gravity model cannot reproduce at the same time 
lattice data of the equation of state of the gluon plasma at very high temperatures, and 
close to the phase transition. This means that fixing the gravity constant G§ from the ideal 
gas limit seems not to be consistent with thermodynamics close to the phase transition. 
In this sense, there is the possibility that the ideal gas limit doesn't correspond to the 
limit of the black hole gravity theory at high temperatures. Is it possible that the gravity 
theory allows more degrees of freedom at high temperatures? Or is the simulation of higher 
terms in the string coupling a' in the gravity action incomplete? Since the agreement of 
the velocity of sound, the spatial string tension and the Polyakov loop in the string frame 
is better, the question arises whether the gravity action approximates the string action 
truthfully. We will further address this problem, and analyze possible solutions [31] in 
forthcoming work. 

The approach presented here and in previous references, see e.g. Refs. JTOj [11], is a 
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phenomenological gravity theory motivated by non-critical string theory. The results are 
subject to 0(1) a' corrections and one can only hope that they capture the expected j3- 
function behavior. It has to be mentioned also as a caveat that the contribution of the 
coupling of the world-sheet to the dilaton field may very well change the quantitative, and 
even the qualitative result substantially. This should be checked in future works. 
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Appendix A: Numerical solution of Einstein equations for the black hole 
metric 

In this appendix we discuss in details the procedure to solve numerically the system of 
Einstein equations given by Eqs. ()3"8l-( l4"Tl) . A numerical solution of the system demands a 
good starting point. The boundary at z = has the disadvantage that b(z) is singular at 
this point. The horizon at z = Zh is not a good expansion point either, since the inverse 
of the black hole factor is singular there. Practically it is possible to start at some value 
close to the horizon, z% = Zh — e- The initial values can then be expanded in terms of their 
values at the horizon as 

1, 



b{zi) 



-W h e + ^e 2 + ... 



h + 7$W h e + 



2 /8 
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:b h W 2 h 



WL)ble 2 + 
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l-y/b^e+- [b h W^--blW h VK + -W^I^ ) e 2 + 
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f'{Zi 



4vrTe + U 2 h V h e 2 + j^T (h h W 2 h + W') b h e 3 + 
= -4ttT - b\V h e - ItiT (h h W 2 + W') b h e 2 + .. 



where we use the notation 

b h = b{z h ) , a h = a{z h ) , W h = W(z h ) , V h = V(a h ) . 



(81) 
(82) 

(83) 

(84) 
(85) 

(86) 



From Eqs. ( 1541 . (I85j) and f(zh) = one can derive Wh and its derivatives. The expression 
for W' h follows from Eq. (1381 by applying l'Hopital rule in the second term of the r.h.s. To 
compute the expression of the second derivative W%, one derives Eq. (|38|) with respect to 



z once, and uses the result of W' h . The result is 
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w' h > 
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where we have defined for simplicity of notation 
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and 

h dz 



• (91) 



The derivative of V with respect to a, i.e. V(a), can be computed analytically from the 
analytical expression of V(a), c.f. Eq. (jSJ). 

Our procedure to solve the system of first order differential equations Eqs. (|38l)-( HT|) 
follows Refs. [T5"| |4~2] . First we choose arbitrary values for the functions at the horizon, 
namely 

&(&) = 10, a(&) = 0.5, (92) 
and the initial values for temperature and e, 

T=lGeV, e=Hr 9 GeV _1 . (93) 

We rewrite the initial conditions, Eqs. ( 18T|) - (!85|) . in the £ coordinate, so that £, = and 
6i — 6 + e = e - The parameter e is chosen very small, such that £, is very close to £/j 
and the initial conditions are accurate enough. The variables £h, 7" differ from 2^, T by a 
rescaling factor. Then one integrates numerically the system to get solutions Wi(£), &i(£), 
ai(£) and /i(£) in some interval £1 < £ < £h, where 61 diverges at £1. Since the system of 
equations fl38l) - fj4Tl) is invariant under three different rescalings [15], one can make use of 
these properties to find a solution which has the right boundary conditions. In step 2, one 
shifts the £ coordinate, so that the ultraviolet divergence of b is at the origin. The new 
solution reads 



W 2 {p) = WMp+£i) , b 2 {p) = 6i(p+£x) , a 2 (p) = «i(p+£i) , h(p) = A(p+£i) . 

(94) 

In step 3, one chooses 6/ in such a way that f^(p = 0) = 1, which is the correct value 
of / in the UV, 



W 3 {p) = W 2 (p)^/5~ f , h{p) = b 2 {p)/y/T f , a 3 {p) = a 2 (p) , f 3 (p) = f 2 (p)/S f . 

(95) 

Finally one rescales 63: 



W(z) = W 3 (zS b ) , b(z) = h(zS b )S b , a(z) = a 3 (z5 b ) , f(z) = f 3 (z5 b ) . (96) 
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The value 5b is determined by comparing the zero temperature metric bo(zo) and the 
finite temperature metric b 3 (z) at a point z which has the same ultraviolet coupling as the 
zero temperature solution. In practice we choose ao(zo) = 0.07. Then the rescaling factor 
5b is given as 

1 = a Oo)/«3(^) ( 97 ) 
5 b = h{z )/h{z) (98) 
z = z 5 b . (99) 
Starting from the values of b(£h), a(£h) an d T mentioned above, one gets 

& = -0.334 GeV^ 1 , 5 f = 0.960 , 5 b = 0.353 . (100) 

By this procedure the QCD-parameter A in the asymptotic logarithms of both solutions 
also agree. Because f'(zh) = —4tiT, the operations performed in Eqs. ( |94l) -( |96|) rescales 
the value of T to the right temperature T, 

t = ^r. (101) 

From Eqs. fllOOp . ( llOip and the value T = 1 GeV one gets T = 368 MeV. In practice we 
solve the equations of motion for different temperatures by choosing different values for 
a(£/i). Note that a(£h) is invariant under the set of rescaling Eqs. (J94l-fl96l). 



Appendix B: Ultraviolet properties of the thermodynamic quantities 

We will study in this appendix the analytical properties of the thermodynamic quanti- 
ties in the ultraviolet. % Being the UV expansion of b{a) given by (c.f. Eq. ( !29l ) 



z(a) 



1 - ^/3o« + |r (22/3 2 - 9/3i) a 2 - -i- (602/3 3 - 540/3 /3i + 81&) a 3 + 0(a 4 ) 
y oi zioi 

(102) 

the entropy can be computed by evaluating the above expression at the horizon, i.e. by 
computing b(ah). The information on the horizon is given by the function f(z), and so 
one should study its behavior in the UV. Combining Eqs. (1391) and (j41j) one gets 

f" b' 

J + 3 T = - (103) 



3 In this analysis we will explicitly show every expression including all the perturbative orders one needs 
to compute the thermodynamic quantities up to 0(a 3 ). 
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Given b(z) one can solve this equation to get f(z). Two integration constants are needed, 
which as usual are chosen by imposing the boundary conditions /(0) = 1 and f(z h ) = [TO] . 
The result for the UV asymptotics of f(z) is 



f(z) 



4 Qj( z h) 



(104) 



where 

QA 



z) = l + ^f3 a(z) - 1(7(3% - 6^)a 2 (z) + 



9 



162 



(271/3 3 - 396/3 /3i + 72/3 2 )a 3 (z) 



0(a 4 ) , 
(105) 

and a(z) is given by Eq. ( 1261) . One can arrive at this result by considering the explicit 
expression of b(z) := b(a(z)) given by Eq. (11021) . insert it into Eq. f |103[) and solve the 
equation reexpressing the result in powers a(z). A much easier way to arrive at this result 
is to consider a general scheme f(z) = 1 — (z A /Ch) ■ (1 + f\ot{z) + / 2 a 2 (z) + f 3 a 3 (z) + ...), 
and then introduce it and Eq. (11021) into Eq. (11031) . The derivative of a(z) is given by 



da 

dz 



(106) 



This useful relation can be proved easily from the explicit expression of a(z), c.f. Eq. 
Using Eq. (I106j) . the equation of motion Eq. (11031) can be expressed in powers of a(z) 
one can easily identify the coeficients Ch, fi, f%, fz, ■ ■ ., which fulfill the equation, 
result is given by Eqs. (11041) and (11051) . 

From Eq. ( 11041) and using Eq. ( 11061) . the derivative of f(z) then evaluates to 



and 
The 



/'(*) 



-Az z 



4Qf( z h) 



1 + |/3 «(^) - |(2/5 ( 2 



+^r(26/3 3 



- W x )a\z) 
36/3 /3i + 9/3 2 )a 3 (z) + C(a 4 ; 



The temperature is obtained by evaluating the above expression at the horizon 



T 



/'(**) 
Air 



KZ h 



1 + 



Pi 



a. 



6 



(5/3 2 - A^)al + 0{a\) 



(107) 



(108) 



We have corrected some error in the computation of Ref. [10J, where the authors consider 
a factor —4/9 instead of 1/3 at order O(o%) in the bracket of Eq. (1108D . 



The entropy density easily follows by evaluating Eq. ( I102j) at the horizon, and using 
the relation between Zh and T given by Eq. (11081) . Then one gets 



s(T) 



4G, 



-b\z h ) 



4G7 



rp3 



1 -^oah + ^(H/3o-6/3i) a 2 h 



3 
1 
162 



9 



(163/3 3 - 252/3 /3i + 72/3 2 ) a\ + 0(a\ 



,(109) 
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which corresponds to the weak coupling expansion of the entropy. To deal with Eq. (I109p 
one needs to know the temperature dependence of an- Taking into account the relation 
between T and Zh given by Eq. (II 08ft . it can be proved thatQ 

a h = a T + ^4 + 0(4) , a T = a(z= , (110) 

where «y is defined as indicated. One can prove that by using the explicit expansion given 
by Eq. (I26I) . By considering z — > l/(jrT) in Eq. (11061) one gets@ 

,d&T,h n 2 n 3 ( 4 Q 3 



T-^ = -p a* T>h - ^a 6 Tih - + ftj a^ h + 0{a^ h ) . (Ill) 

The ah version of this formula easily follows by considering Eq. (IllOp . Eq. (Illip is very 
useful, and it can be used for instance to compute easily the pressure from Eq. dl09[) (see 
Section [6] for a discussion). The energy density follows trivially from Eqs. (11091) and ( 1621) 

e(T) s p 3n 3 £ 3 



T 4 T 3 T 4 16G, 



l--P o a h +-(20Z-p 1 )aj i 



(112) 



Note that in the weak coupling expansion of the thermodynamics quantities, Eqs. ( 11 09j) . 
( 162|) and f 1 1 1 2 [) . there are no half-integer powers in a, i.e. a 3//2 , a 5//2 , . . . , as we don't consider 
loops contributions, in contrast to the weak coupling expansion in QCD [3D]. These results 
extend to O(a^) the results of Ref. [TDJ. From the energy density and pressure one can 
compute the trace anomaly. It reads 



A(T) 



(3(a) (TrF 2 ,) e - 3p 



7T 



i£ 3 



4a 2 T 4 



y4 



12G, 



/3 2 a 2 - -A, (2/3 2 - 3/30 a* + ( Q 



113) 



The trace anomaly is related to the vacuum expectation value of the gluon condensate. As 
it was pointed out in Ref. [10] and discussed by us in Sec. 131 the gluon condensate appears in 
the UV expansion of the difference between the black-hole and zero temperature solutions, 
c.f. Eqs. P6|) - fj4*9|) . In this Appendix we have not taken into account power corrections 
in z. However, just by computing b 3 (zh) using Eq. fjlBT) . it is straightforward to prove 



4 Note that because of Eq. (|110[) the expressions of the thermodynamics quantities up to 0(a h ) remain 
valid when one substitutes ah by ay. 

5 The definition of ax given in Eq. (|110[) doesn't agree with the usual definition of the running coupling 
at finite temperature, for which the prescription E ~ irT is usually taken. Both prescriptions differ at 



O(a^), i.e. cce — ar + 0(a\,), as it is evident from Eq. (|1 1 1 1) . 
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that the correction ~ z A induces a contribution ~ 1/T 4 in s(T)/T 3 , p(T)/T 4 , e(T)/T 4 and 
A(T), and so this corresponds to much a lower order contribution in the UV expansion 
performed previously. By the same way, the power correction in a(z) induces a correction 
~ 1/T 4 in a(zh) which is subleading in our UV analysis, and so it is enough to identify 
a(z h ) with a (zh) at this level, as it has been done in previous formulas. 

From previous analysis we can derive easily the expressions for the weak coupling 
expansion of the specific heat per unit volume 



T 



d 2 p 3tt 

W 1 



3/>3 



4G E 



rp6 



1 - ~/3 a h + ~ (5/3 2 - 2fa) a 



162 



(-295/3^ + 396/3 /3i - 72/3 2 ) a\ + 0(a 



;ii4) 



and speed of sound 



s 

( v 



\&«h + l& W - 4/?i) a 3 h + 0(ai) 



(115) 



Eq. flTTD follows by using Eqs. f lTTT]) and while Eq. f TTTol) is obtained from Eqs. <HJJU\i 
and fTTTD . 

For completeness of this apendix, we study the high temperature behavior of the 
Polyakov loop. At the end one wants to express the result as an expansion in powers 
of ah, and the easiest way to proceed is to work in coordinates dependent on the running 
coupling a as a variable, instead of z. The relation between z and a is given by [1] 



da 
dz 



-b(a)e 



-D(a) 



;ii6) 



where the function e D ^ reads 



(3(a 



exp 



3a 2 



da 



117) 



Using Eqs. ( 1116)) . one can rewrite the Nambu-Goto action which involves an integration 
in z, Eq. ( |75l) . as 



Sk 



H dae D{a) b(a)J . 



'118) 



2nP s T 

In the intermediate steps we will make use explicitly of the UV /3-function up to 4-loops 
order just for completeness, i.e. 



(3(a) = -(3 a 2 - /3i« 3 - (3 2 a 4 - (3 3 a 5 H , 



(119) 
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but our final result of the Polyakov loop up to 0(a^^ 3 ) will depend only on ft and ft, c.f. 
Eq. (fI23|) . Inserting the UV ft function, Eq. (HUD, into Eq. ffTTTj) . one gets 



3 D(«) 



fta 2 \9 ft/ « 81 V A) ft 



(120) 



Note that Eq. (11161) combined with the expansions of 6(a) and e D ^ given by Eqs. ( 11021) 
and (11201) respectively, leads to Eq. ( I106p . The main difficulty is to express b as a function 
of a. The UV expansion of 6(a) is given by Eq. (I102p . To compute z(a) one has to invert 
Eq. (T26D. The result is 



1 



A(fta)/V/3 2 



exp 



fta 



4 



ft 2 fa 



+ (^o 2 - §A - 

From Eqs. (1TP2]) and ( ETJ , one gets 



2ft 4 



(ft 3 - 2ftftft + /3o/5s) « 2 + 0(a 



'121) 



6(a) 



pi , 



1 + ^^ a 



+ — g (ft 4 - 2ftftft(ft 2 + ft) + ft 2 (ft 3 + ft 2 + ft 2 ft)) a 2 + 0(a 3 



(122) 



Then inserting Eqs. fl 1 2 [) and fl 1 22[) into Eq. (11181) . and performing the integration in a, 
one arrives at 



- ^4 (l + ^fa*h + ^ (161ft 2 + 72ft) a 2 + 0(a 3 h/ 



(123) 



where S^ G is divergent and comes from the lower limit in the integration. To arrive at 
Eq. ( 11231) one has to make use of Eq. ( I12ip evaluated at the horizon, and use the relation 
between T and Zh given by Eq. ( 11081) . Then the renormalized vacuum expectation value 
of the Polyakov loop writes 



Lr(T) 



e NG = exp 



r 02 



_«3 l 1+ _fc ah + _ ( 161/3 2 + 72/3 J ^ + 



(124) 



where S™ G = S NG — S® G . Note that Lr(T) tends to 1 from above in the high temperature 
limit, which is the behavior predicted by standard perturbative QCD. 
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